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Abstract 

A  dynamic  model  of  a  centrifugal  compressor  capable  of  system  simulation  in  the  virtual  test  bed  (VTB)  computational  environment  is  presented. 
The  model  is  based  on  first  principles,  i.e.  the  dynamic  performance  including  the  losses  is  determined  from  the  compressor  geometry  and  not  from 
the  experimentally  determined  characteristic  performance  curves.  In  this  study,  the  compressor  losses,  such  as  incidence  and  friction  losses,  etc., 
are  mathematically  modeled  for  developing  compressor  characteristics.  For  easy  implementation  in  the  VTB  platform,  the  non-linear  governing 
equations  are  discretized  in  resistive  companion  (RC)  form.  The  developed  simulation  model  can  be  applied  to  virtually  any  centrifugal  compressor. 
By  interfacing  with  a  composite  system,  such  as  a  Brayton  cycle  gas  turbine,  or  a  fuel  cell,  the  compressor  dynamic  performance  can  be  evaluated. 
The  surge  line  for  the  compressor  can  also  be  determined  from  the  simulation  results.  Furthermore,  the  model  presented  here  provides  a  valuable 
tool  for  evaluating  the  system  performance  as  a  function  of  various  operating  parameters. 
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1.  Introduction 

Analytical  performance  prediction  method  plays  an  impor¬ 
tant  role  in  designing  a  centrifugal  compressor  by  way  of  pre¬ 
dicting  the  overall  dimensions  and  performance  curve  of  the 
stage.  First  of  all,  the  analytical  method  can  be  used  to  perform 
parametric  studies  to  demonstrate  the  influence  of  changes  in 
geometry  on  the  performance  under  both  design  and  off-design 
conditions.  Therefore,  the  availability  of  reliable  analytic  tool 
saves  expensive  experimental  development.  In  addition,  analytic 
models  for  predicting  the  overall  performance  of  the  compressor 
can  be  effectively  used  in  predicting  the  overall  performance  of 
a  combination  system  where  the  compressor  is  a  component  of 
such  a  system. 

In  order  to  derive  the  analytical  tools,  the  losses  occurring 
throughout  the  stage  must  be  specified.  The  origin  and  effects 
of  loss  mechanism  were  discussed  in  details  in  references  [4-6]. 
These  papers  presented  prediction  methods  based  on  modeling 
various  losses  throughout  the  stage.  According  to  references 
[1,2,7],  the  losses  in  the  centrifugal  compressor  are  classified 
into  incidence  loss,  friction  loss,  clearance  loss,  backward  loss 
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and  volute  loss.  Reference  [3]  presents  comparative  study  of  the 
existing  modeling  techniques  and  their  limitations  for  predict¬ 
ing  incidence  loss.  There  are  two  widely  used  models:  one  is  the 
constant  pressure  incidence  model  and  the  other  is  the  NASA 
shock  loss  theory  as  described  in  references  [2,3].  In  the  current 
study,  the  NASA  shock  loss  theory  is  used.  It  can  be  mentioned 
here  that  very  little  difference  exists  between  the  two  models  for 
centrifugal  compressor  as  expounded  in  reference  [1].  As  to  the 
friction  loss,  references  [2,3]  apply  the  energy  and  momentum 
equations  to  pipe  flow  with  surface  frictions  to  get  the  loss  coef¬ 
ficient  related  with  Reynolds  number.  For  this  study,  we  used 
the  approach  of  references  [2,3]  to  predict  the  friction  loss. 

The  virtual  test  bed  (VTB)  provides  an  effective  computa¬ 
tional  environment  to  simulate  the  dynamic  performance  of  the 
centrifugal  compressor  [12].  The  non-linear  model  equations 
based  on  energy  transfer  are  discretized  in  resistive  compan¬ 
ion  (RC)  form  for  effective  implementation  in  the  VTB  plat¬ 
form.  After  the  compressor  model  is  developed  and  validated;  a 
dynamic  simulations  of  a  fuel-cell  system  assembled  in  the  VTB 
environment  are  carried  out  to  analyze  the  compressor  perfor¬ 
mance  along  with  the  system  performance.  Details  of  VTB  is 
given  in  reference  [11,13]. 

The  remainder  of  this  paper  is  organized  as  follows.  The 
mathematical  identify  equations  for  the  centrifugal  compressor 
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Nomenclature 

a 

hydraulic  perimeter  (m) 

Ai 

absolute  inlet  velocity  (ms-1) 

cp 

specific  heat  (kgkJ-1  K-1) 

Cn 

inlet  radial  velocity  of  fluid  (ms-1) 

Cr2 

outlet  radial  velocity  of  fluid  (ms-1) 

Cl 

absolute  inlet  velocity  (ms-1) 

C2i 

the  loss  of  outlet  tangential  velocity 
(ms-1) 

of  fluid 

Cei 

inlet  tangential  velocity  of  fluid  (ms-1) 

ce2 

outlet  tangential  velocity  of  fluid  (ms-1 

) 

D 

mean  hydraulic  channel  diameter  (m) 

Aii 

hub  casing  diameter  (m) 

Ai 

inducer  tip  diameter  (m) 

D\ 

inlet  average  diameter  (m) 

d2 

outlet  average  diameter  (m) 

f 

friction  factor 

m 

mass  flow  rate  (kg  s-1) 

J 

spool  moment  of  the  inertia  (kg  m2) 

l 

mean  channel  length  (m) 

PO 

inlet  pressure  of  compressor  (Pa) 

P2 

outlet  pressure  of  compressor  (Pa) 

n 

inlet  average  radius  (m) 

r2 

outlet  average  radius  (m) 

Re 

mean  Reynolds  number 

To 

inlet  temperature  of  compressor  (K) 

t2 

outlet  temperature  of  compressor  (K) 

u 

blade  tangential  velocity  (ms-1) 

Ui 

the  inlet  tangential  velocity  of  the 
(ms-1) 

impeller 

u2 

the  outlet  tangential  velocity  of  the 
(ms-1) 

impeller 

Wc 

power  delivered  to  the  fluid  (J) 

W\ 

ideal  inlet  relative  velocity  (ms-1) 

Wih 

actual  inlet  relative  velocity  (ms-1) 

Wei 

the  loss  of  relative  inlet  tangential 
(ms-1) 

velocity 

OL2 

fluid  inlet  angle  (°) 

^2b 

blade  outlet  angle  (°) 

/?i 

fluid  inlet  angle  (°) 

Alb 

blade  inlet  angle  (°) 

A/ibf 

backflow  loss  (Jkg-1) 

A/ici 

clearance  loss  (Jkg-1) 

A/idf 

friction  loss  for  diffuser  (Jkg-1) 

Ahf 

friction  loss  (Jkg-1) 

Ah[ 

incidence  loss  (Jkg-1) 

A/^id 

incidence  loss  for  diffuser  (Jkg-1) 

A/?ideal 

ideal  specific  enthalpy  delivered  to 
(Jkg”1) 

the  fluid 

AJl[f 

friction  loss  for  inducer  (Jkg-1) 

Aha 

incidence  loss  for  inducer  (Jkg-1) 

Ahy 

volute  loss  (Jkg-1) 

A??bf 

efficiency  loss  increases  with  backflow 

A)?ci 

efficiency  loss  increases  with  clearance 

A)7v 

efficiency  loss  increases  with  volute 

y  specific  heat  ratio 

r/i  efficiency 

poi  constant  stagnation  inlet  density  (kg  m-3) 
a  slip  factor 

rc  compressor  torque  (N  m) 

rt  drive  torque  (N  m) 

co  rotational  speed  (rpm) 


are  presented  in  Section  2.  The  RC  model  formulation  is  derived 
in  Section  3.  Section  4  gives  the  simulation  results  of  compres¬ 
sor.  In  Section  5,  discussion  of  the  performance  and  analysis 
of  the  compressor  in  the  overall  fuel-cell  system  are  presented. 
Conclusions  are  made  in  Section  6. 


2.  Identify  equations 


Energy  transfer  is  mainly  considered  in  developing  compres¬ 
sor  characteristic.  The  fluid  enters  the  compressor  rotor  at  one 
radius  with  one  velocity  and  leaves  at  another  radius  with  another 
velocity.  The  change  in  momentum  of  the  fluid  is  derived  from 
the  work  done  by  the  rotating  rotor,  which  is  driven  by  an  exter¬ 
nally  applied  torque.  The  compressor  characteristic  equation  can 
be  described  as  follows: 

1 

=  -  •  (T t  -  Tc)  (1) 


CO  ■  Tc  =  Wc  =  rn  ■  A/iideal 


P  2 


mgih  U)  •  Ah  ideal 

To  ■  Cp 


y/(y- 1) 

•  Po 


Wc  o 

A /l ideal  =  —  =  O  ■  lj\  (4) 

m 

Eq.  (1)  presents  momentum  torque  equation.  The  difference 
between  external  drive  torque  rt  and  compressor  torque  rc  used 
to  change  the  fluid  momentum  is  expended  to  cause  the  rotor  to 
rotate. 

Eq.  (2)  comes  from  the  Euler’s  pump  equation,  where  m  is  the 
mass  flow  rate  of  the  compressor.  The  specific  enthalpy  A/zideai 
delivered  to  the  fluid  should  be  equal  to  the  power  delivered  to 
the  fluid  Wc. 

Eq.  (3)  is  derived  from  isentropic  compression  relations 
between  pressure  and  temperature.  Various  kinds  of  losses  are 
considered  in  determining  the  efficiency  iji  U,  which  is  a 
function  of  the  mass  flow  rate  and  rotational  speed). 

Eq.  (4)  further  specifies  the  specific  enthalpy  A/zideai,  where 
we  assume  that  there  are  no  stationary  pre-whirl  vanes  and  air 
approaching  the  impeller  does  not  have  a  tangential  component 
of  velocity,  a  is  called  the  slip  factor,  the  ratio  between  out¬ 
let  tangential  velocity  of  the  fluid  Cq 2  and  the  outlet  tangential 
velocity  of  the  impeller  C/2.  For  the  details  of  the  derivation  of 
Eq.  (4),  please  refer  to  Appendix  A. 
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2.1.  Loss  modeling 


2.1.1.  Incidence  losses 

Incidence  losses  come  from  the  off-design  velocity  triangles 
at  the  impeller  eye  (inducer)  causing  flow  separation  and  are 
therefore  at  the  design  point  this  loss  is  zero. When  the  fluid  is 
entering  the  impeller  eye  at  a  relative  velocity  of  W\  and  with 
an  inlet  angle  p i,  as  shown  in  Fig.  1,  it  is  assumed  that  the 
fluid  instantaneously  changes  its  direction  to  comply  with  the 
blade  inlet  angle  pi^.  As  a  result,  the  relative  velocity  of  the 
fluid  is  changed  from  W\  to  Wib.  According  to  references  [2,3], 
NASA  theory  it  is  assumed  that  the  kinetic  energy  associated 
with  tangential  component  Wei  is  destroyed  as  the  fluid  adapts 
to  the  blade  direction.  Thus,  the  energy  loss  due  to  incidence  is 
given  by 


From  Fig.  1,  it  can  be  shown  that 


cos  — 


U l  ~  Cei 
W\ 


and  sin  p\ 


Crl 

Wi 


where  Cei  and  Cv\  are  the  inlet  tangential  velocity  and  inlet 
radial  velocity  of  fluid,  respectively,  and 


sin(/3ib  —  P\) 

Wei  =  V*  PiVi  =  (cos^i  -  cotytfib  sin  pi) •  Wi  (7) 
sm/3ib 


From  Eqs.  (6)  and  (7) 


Wei  =  U i  —  Cei  -  cot^ib  cr  1 


So  the  incidence  loss  can  be  written  as 

1  9  1  /  cot Pib  m\2 

Aha  =  -(Ui  -  Cei  -  cot/Jlb  crl)2  =  -  £/i - 

2  2\  poiAi  J 

(9) 


Similar  to  the  flow  situation  in  the  impeller,  when  the  fluid  enters 
the  vane  diffuser  at  a  velocity  of  C2  and  with  an  inlet  angle  012, 
as  shown  in  Fig.  2,  it  is  assumed  that  the  fluid  instantaneously 
changes  its  direction  to  comply  with  the  diffuser  inlet  angle 
o'2b-  As  a  result,  the  velocity  of  the  fluid  is  changed  from  C2  to 
C2b-  According  to  reference  [1,10],  the  kinetic  energy  associated 


with  the  tangential  component  C2i  of  the  velocity  is  lost.  So  the 
incidence  loss  can  be  expressed  as 

C 2 

Aftid  =  -f  (10) 

From  Fig.  2,  we  can  obtain 

1  2 

A /lid  =  ~(crU2  ~  cot  a  2b  Cr2)  (11) 

where  Cr2  is  the  outlet  radial  velocity  of  the  fluid  and  C02  is  the 
outlet  tangential  velocity  of  fluid. 

For  simplicity,  we  further  assume  that  Cr  1  =  Cr2. 

U\  =  Cri  cotjftib  =>  Cr2  =  U\  tan  /3ib  (12) 


From  Fig.  2,  it  follows  that 
Cr2  U\  tan  p lb 


tan  o'2b  — 


So 


Ce  2 


oU2 


«2b  =  a  tan 


Of  tan  Pi\y 

0D2 


(13) 


(14) 


and  according  to  the  above  equations  the  diffuser  incidence  loss 
can  be  written  as 


A/zid 


1  foD2Ui 

2  \  D\ 


m  •  cot  a 2b  \ 
P10M  J 


(15) 


2.1.2.  Friction  losses 

According  to  reference  [2],  the  friction  loss  in  the  impeller 
can  be  defined  as 

/  W?u 

A%  =  4/-^  (16) 

where  /  is  the  friction  factor  for  appropriate  mean  Reynolds 
number,  /  the  mean  channel  length  and  D  is  the  mean  hydraulic 
channel  diameter,  which  is  defined  as 

/  =  0.3164(/?e)-0,25  (17) 


Fig.  2.  Incidence  losses  of  fluid  entering  diffuser. 
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The  mean  hydraulic  channel  D  is  defined  as 


a 

From  Fig.  2,  it  is  seen  that 

Wib  =  Wi 

sin/Ji  sin /Jib 

and  using  sin  (i\  =  (Cr\  /  W\ ),  we  can  get 


(IB) 


(19) 


Wib  = 


C  i 


sin  /Jib 

Finally,  we  can  get  the  friction  loss  at  impeller 


(20) 


Ah\f  =  kfjh  = 


Chl 


2 Dp\A\  sinz^ib 


m 2 


(21) 


For  diffuse  the  friction  loss  can  be  modeled  in  a  similar  manner 
as  in  the  impeller 


A/zdf  =  kf&m 


2 


chl 

2 Dp\A\  sin2  /3ib 


(22) 


2.1.3.  Other  losses 

According  to  reference  [2],  the  clearance  loss  is  a  function 
of  the  clearance  to  passage  width  at  the  tip.  Pampreen  showed 
that  the  stage  efficiency  loss  increases  with  clearance  and  can  be 
expressed  as  A^d  ~  0.3(/ci/£>2),  where  lc\  is  the  tip  axial  clear¬ 
ance  and  Z?2  is  the  impeller  tip  width.  For  back  flow  loss,  no 
theory  or  mathematical  model  exists  at  present  to  describe  the 
backflow  loss,  Dean  and  Young  suggested  a  backflow  loss  of 
three  points  of  stage  efficiency  as  typical:  A^bf  =  0.03.  Accord¬ 
ing  to  reference  [1],  the  volute  loss  is  assumed  to  lie  within 
two  to  five  point  of  efficiency:  0.02  <  Arjv  <  0.05.  According  to 
simulation  results  of  this  paper,  under  the  operating  conditions 
of  a  rotational  speed  of  35,000  rpm  and  in  the  neighborhood  of 
the  surge  line,  the  relative  magnitudes  of  clearance,  backflow 
and  volute  losses  are  13.33,  12.5  and  14.58%,  respectively,  of 
the  total  loss.  It  should  be  noted  that  the  relative  magnitudes  of 
clearance,  backflow  and  volute  losses  will  decrease  gradually  as 
mass  flow  rate  goes  up,  this  is  due  to  the  fact  that  at  higher  mass 
flow  rates  the  incidence  loss  and  friction  loss  will  increase. 


2.2.  Efficiency 


The  isentropic  efficiency  of  the  compressor  is  defined  as 

A  h  ideal 


U)  = 


A/Zideal  T  Ah\ 


-  Arid  ~  A?7bf  -  A77, 


(23) 


OSS 


where  in  this  paper,  the  loss  Ah\oss  is  equal  to  the  sum  of  inci¬ 
dence  losses  and  friction  losses. 


A/^ioss  —  A/iji  +  Ah(^[  +  Ah[f  -I-  A/idf 


(24) 


3.  RC  model  formulation 


The  resistive-companion  method  [8,9]  provides  a  way  to 
account  for  natural  conservation  laws  by  defining  a  pair  of 
across  and  through  variables  at  each  terminal.  The  device  object 


interacts  with  the  VTB  network  solver  by  providing  the  device 
conductance  matrix  and  the  history  vector  at  each  simulation 
time  steps,  so  that  the  solution  to  the  entire  circuit  can  be  sought 
by  the  solver  based  on  the  mathematically  equivalent  nodal  cir¬ 
cuit  analysis.  The  solver  requires  that  the  relations  of  the  terminal 
variables  for  each  device  be  written  in  the  following  standard 
form: 

7(0  =  G  •  V(t)  -  Bit  -  h )  (25) 

where  I(t)  is  the  through  variable  vector,  V(t)  the  across  variable 
vector,  G  the  conductance  matrix,  B(t  —  h )  the  history  vector  of 
the  device  and  h  is  the  simulation  time  step.  Notice  that  although 
the  term  “conductance  matrix”,  inherited  from  electric  network 
analysis,  is  used,  the  terminal  variables  are  more  generally  across 
and  through  variables,  not  necessarily  voltage  and  current. 

In  Fig.  3,  it  is  shown  that  the  compressor  interacts  with  the 
external  world  through  its  two  terminals.  Node  0  is  a  mechanical 
terminal  where  torque  r  is  through  variable  and  rotational  speed 
co  is  an  across  variable,  through  which  compressor  connects  with 
outer  driving  machine  such  as  a  motor  or  a  gas  turbine.  Node  1 
is  a  fluid  terminal  where  mass  flow  rate  m  is  a  through  variable 
and  pressure  p  is  an  across  variable,  by  which  compressed  fluid 
is  outputted. 

Eq.  (1),  co  =  j  •  (rt  —  rc),  is  discretized  within  one  time  step 
using  trapezoidal  method. 

2  J  2  J 

Tt (0  =  —00{t)  +  Tc(t) - —Co(t  ~  h ) 

h  h 

-  Tft  -  h)  +  x c(t  -  h)  (26) 

For  the  non-linear  equations  (2)  and  (3),  Newton-Raphson  iter¬ 
ations  is  introduced  to  solve  these  equations.  By  eliminating  the 
intermediate  variables  rc,  771  (m,  U)  and  A/zideal,  we  can  obtain 
standard  RC  equations  presented  below: 


Fig.  3.  Centrifugal  compressor  icon  in  VTB. 
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where  the  history  for  m  is 

dm  dm 

bo(t  —  h)  =  —m(t  —  h)  +  - —  •  pi(t  —  h)  +  •  co(t  —  h) 


dp2 


dco 


(28) 


3rt  3rt 

b\ (t  —  h)  =  —T\(t  —  h)  +  - —  •  p2(t  —  h)  +  —  •  co(t  —  h ) 


d/?2 


(29) 


Finally,  we  can  obtain  the  complete  RC  model  equations 
(27)-(29)  for  the  centrifugal  compressor. 


4.  VTB  simulation 


An  integrated  hybrid  solid  oxide  fuel  cell  energy  system, 
assembled  in  the  VTB  environment,  is  shown  in  Fig.  4.  The  cen¬ 
trifugal  compressor  is  driven  by  a  motor,  and  the  compressed  air 
is  supplied  to  fuel  cell.  The  system  consists  of  a  fuel  compressor, 
an  air  compressor,  a  fuel  cell,  a  combustor,  a  gas  turbine,  a  heat 
exchanger  and  a  fuel  tank.  The  compressed  air  is  pre-heated  by  a 
heat  exchanger  prior  to  entering  the  fuel  cell,  where  the  natural 
gas  is  internally  reformed.  After  the  electrochemical  reaction, 


electrical  energy  is  produced  which  is  associated  with  simulta¬ 
neous  heat  generation  during  the  process.  For  further  extraction 
of  energy  the  un-reacted  high-temperature  gases  are  channeled 
to  the  combustor  for  complete  combustion  there.  The  exhaust 
gas  then  expands  through  the  gas  turbine  whereby  mechanical 
power  is  generated.  The  exhaust  gas  discharging  from  the  gas  tur¬ 
bine  is  further  used  for  the  purpose  of  regeneration,  by  passing  it 
through  the  heat  exchanger  where  it  pre-heats  the  compressed  air 
going  to  the  fuel  cell  stacks.  The  load  of  the  compressor  depends 
on  the  number  of  cells.  The  output  of  this  fuel  cell  system,  in 
the  form  of  electrical  power,  is  supplied  to  an  electrical  circuit 
load.  It  can  be  mentioned  here  that  for  the  fuel-cell  system  the 
compressor  had  a  flow  rate  of  less  than  1  kg  s_1,  pressure  ratio 
of  less  than  4,  and  the  maximum  compressor  power  of  332  kW 
and  power  generated  by  the  overall  system  was  less  than  1  MW. 
The  modeling  details  for  the  fuel  cell  and  other  components  are 
not  presented  in  this  paper  because  the  focus  of  this  paper  is  the 
centrifugal  compressor  modeling. 

The  physical  dimensions  of  some  of  the  parameters  of  the 
centrifugal  compressor  used  in  this  system  are  listed  in  Table  1 . 

The  simulation  results  obtained  in  VTB  are  presented  as  four 
representative  curves,  as  shown  in  Fig.  5,  including  the  mass  flow 


Table  1 

Parameters  for  dynamic  centrifugal  compressor  model 

Fig.  4.  Simulation  system. 

Parameter 

Value 

Parameter 

Value 

Slip  coefficient 

0.9 

Outlet  tip  diameter  (m) 

0.15 

Blade  inlet  angle  (°) 

45 

Outlet  hub  diameter  (m) 

0.10 

Diffuser  inlet  angle  (°) 

45 

Mean  impeller  channel  length  (m) 

0.12 

Inlet  tip  diameter  (m) 

0.07 

Mean  diffuser  channel  length  (m) 

0.1 

Inlet  hub  diameter  (m) 

0.034 

Spool  moment  of  inertia  (kg  m2) 

5 

Pressure  atmosphere  (Pa) 

101,000 

Temperature  atmosphere  (K) 

298 

Density  of  the  fluid  (kg  m~3) 

1.205 

Kinetic  viscosity  of  the  fluid  (m2  s-1) 

0.001 
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Fig.  5.  Dynamic  simulation  results  obtained  in  VTB. 


rate,  output  pressure,  rotational  speed  and  efficiency,  each  of 
which  is  displayed,  respectively,  as  a  function  of  time  in  seconds. 
The  rotational  speed  and  efficiency  curves  start  from  zero  and 
reach  a  steady  state  at  the  end.  Meanwhile,  the  output  pressure 
increases  from  the  atmospheric  pressure  to  a  steady-state  value. 
The  magnitude  of  mass  flow  rate  is  mainly  decided  by  the  fuel 
cell. 

5.  Discussion 

5.7.  Compressor  characteristics 

Further  derivation  based  on  the  VTB  simulation  results  leads 
to  compressor  characteristics  curves  with  a  specific  set  of  param¬ 
eters.  The  extensive  characteristics  curves  studies  can  be  per¬ 
formed  by  intentionally  changing  the  compressor  parameters  as 
needed.  The  resulting  analysis  lays  a  solid  foundation  for  the 
behavior  prediction  of  the  compressor. 

Fig.  6a  shows  efficiency  versus  mass  flow  rate  characteristics 
obtained  from  simulation  under  different  rotational  speeds  V, 
ranging  from  20,000  to  50,000  rpm.  The  observations  from  the 
figure  are  listed  as  follows: 

•  The  optimum  mass  flow  rate,  i.e.  best  efficiency  of  the  com¬ 
pressor  changes  a  function  of  rotational  speed.  For  example, 
the  optimum  mass  flow  rate  is  around  0.15  kg  s-1  in  the  case 
of  20,000  rpm  rotational  speed  and  it  changes  to  0.3  kg  s-1  in 
the  case  of  50,000  rpm  rotational  speed. 

•  The  optimum  mass  flow  rate  increases  as  the  rotational  speed 
goes  up. 

•  Under  a  specific  rotational  speed,  in  case  that  the  mass  flow 
rate  is  beyond  the  optimum  value,  the  efficiency  of  the  com¬ 
pressor  begins  to  decrease  as  the  mass  flow  rate  increases. 
And  this  phenomenon  gets  more  pronounced  at  lower  rota¬ 
tional  speeds. 


•  Thus,  based  on  the  first  two  observations,  for  a  given  mass 
flow  rate  in  the  real  case,  the  optimum  characteristics  of  the 
compressor  can  be  achieved  by  adjusting  the  rotational  speed 
N. 

•  Based  on  the  third  observation,  it  is  apparent  that  operating 
under  higher  rotational  speed  leads  to  a  better  efficiency  espe¬ 
cially  for  larger  mass  flow  rates. 

Fig.  6b  presents  outlet  pressure  versus  mass  flow  rate  char¬ 
acteristics  obtained  from  simulation  under  different  rotational 
speeds  V,  ranging  from  20,000  to  50,000  rpm.  The  observations 
from  the  figure  are  listed  as  follows: 

•  The  outlet  pressure  is  a  function  of  the  mass  flow  rate  and 
rotational  speed  N.  Therefore,  the  desired  outlet  pressure  can 
be  achieved  by  regulating  either  the  mass  flow  rate  or  the 
rotational  speed  or  both  of  them. 

•  The  rotational  speed  imposes  a  much  bigger  effect  on  the 
outlet  pressure  than  that  of  the  mass  flow  rate.  As  a  result, 
adjustment  of  the  rotational  speed  works  more  favorably  for 
the  macro-adjustment  of  the  outlet  pressure  while  the  reg¬ 
ulation  of  the  mass  flow  rate  would  be  preferable  for  the 
micro-adjustment  of  the  outlet  pressure. 

Fig.  6c  illustrates  incidence  loss  versus  mass  flow  rate  char¬ 
acteristics  obtained  from  simulation  under  different  rotational 
speeds  V,  ranging  from  20,000  to  50,000  rpm.  The  observations 
from  the  figure  are  listed  as  follows: 

•  Incidence  loss  decreases  at  first  till  it  reaches  minimum 
value,  i.e.  the  on-design  working  point  without  incidence 
loss.  The  existence  of  the  zero  incidence  loss  is  explained 
by  the  fact  that  the  change  of  the  mass  flow  rate  affects 
the  fluid  velocity;  this  change  further  influences  the 
fluid  inlet  angle.  When  the  fluid  inlet  angle  coincides 
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(b)  Mass  flow  rate  (kg/s) 


(d)  Mass  flow  rate  fka/s) 


Fig.  6.  (a)  m  vs.  efficiency  for  various  rotational  speed;  (b)  m  vs.  pressure  for  various  rotational  speed;  (c)  m  vs.  incidence  loss  for  various  rotational  speed;  (d)  m 
vs.  friction  loss  for  various  rotational  speed. 


with  the  blade  inlet  angle,  the  incidence  loss  becomes 
zero. 

•  Different  rotational  speed  leads  to  different  zero  incidence 
loss  point.  Since  the  fluid  inlet  angle  is  a  function  of  both 
mass  flow  rate  and  rotational  speed,  the  mass  flow  rate  needs 
to  be  regulated  accordingly  as  the  rotational  speed  changes 
so  as  to  maintain  the  coincidence  of  the  fluid  inlet  angel  with 
the  blade  inlet  angle. 

Fig.  6d  depicts  friction  loss  versus  mass  flow  rate  character¬ 
istics  obtained  from  simulation  under  various  rotational  speeds 
N,  ranging  from  20,000  to  50,000  rpm.  The  overlapping  of  all 
the  curves  indicates  that  frictional  loss  is  a  function  of  mass 
flow  rate  only  and  that  the  change  of  the  rotational  speed  has  no 
bearing  on  the  friction  losses. 


5.2.  Parameters  analysis 

Parametric  analysis  plays  a  very  important  role  in  the  design 
of  compressors.  With  the  analytical  method  applied  to  paramet¬ 
ric  studies,  it  is  easy  to  demonstrate  the  influence  of  changes  in 
geometry  on  the  performance  under  both  design  and  off-design 
conditions,  thus  providing  the  guideline  for  compressor  opti¬ 
mization  without  expensive  and  time-consuming  experimental 
efforts. 

Fig.  7  demonstrates  outlet  pressure  versus  mass  flow 
rate  (a)  and  efficiency  versus  mass  flow  rate  characteris¬ 
tics  curve  (b)  with  different  slip  factors  while  all  the  other 
parameters  were  kept  unchanged.  As  shown  in  the  fig¬ 
ures,  both  the  outlet  pressure  and  efficiency  decrease  when 
slip  factor  changes  from  0.9  to  0.8.  This  implies  that, 
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Fig.  7.  (a)  m  vs.  pressure  for  two  different  slip  factors  and  (b)  m  vs.  efficiency  for  two  different  slip  factors. 


a  larger  slip  factor  translates  into  better  compressor  effi¬ 
ciency. 

Fig.  8  shows  outlet  pressure  versus  mass  flow  rate  (a)  and  effi¬ 
ciency  versus  mass  flow  rate  (b)  characteristics  curve  at  different 
inlet  and  outlet  diameters.  The  exact  parameters  for  various 
cases  are  listed  in  Table  2.  As  shown  in  Fig.  8a,  with  all  the 
other  parameters  remaining  unchanged,  smaller  diameter  pro¬ 
duces  smaller  outlet  pressure.  As  observed  from  (b),  the  effect 
of  diameter  size  on  the  efficiency  is  weak  when  mass  flow  rate 
is  around  0.2  kg  s-1,  whereas  the  same  becomes  pronounced  as 
the  mass  flow  rate  increases.  In  conclusion,  the  compressor  with 
smaller  diameter  should  not  be  operating  under  off-design  con- 


Table  2 

Parameters  for  cases  1  and  2 


Inlet  tip 

Inlet  hub 

Outlet  tip 

Outlet  hub 

diameter  (m) 

diameter  (m) 

diameter  (m) 

diameter  (m) 

Case  1 

0.07 

0.034 

0.15 

0.1 

Case  2 

0.06 

0.024 

0.14 

0.09 

x  105 


ditions,  where  the  mass  flow  rate  is  expected  to  vary  significantly 
during  the  operation  of  the  compressor,  leading  to  a  degrading 
efficiency. 

Additional  studies  are  performed  to  investigate  the  effect  of 
blade  inlet  setting  angle  on  pressure,  efficiency  and  losses,  as 
shown  in  Figs.  9a  and  b  and  10,  respectively.  Both  the  efficiency 
and  the  pressure  are  lower  for  smaller  setting  angles,  obviously 
the  losses  shown  in  Fig.  10  has  an  opposite  trend.  Given  the 
same  mass  flow  rate  under  off-design  conditions,  the  efficiency 
of  the  compressor  gets  poorer  for  the  smaller  blade  inlet  setting 
angle. 

5.3.  Validation  of  model 

The  proposed  model  was  validated  against  experimental  and 
analytical  results  from  the  literature.  For  experimental  compar¬ 
isons,  the  parameters  for  the  centrifugal  compressor  is  as  shown 
in  Table  3. 

Fig.  1 1  shows  the  comparison  between  the  VTB  simulation 
results  at  three  different  rotational  speeds  with  the  experiment 


Fig.  8.  (a)  m  vs.  pressure  for  different  inlet  and  outlet  diameters  and  (b)  m  vs.  efficiency  for  different  inlet  and  outlet  diameters. 


Pressure  Ratio  Losses  (J/kg) 
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Fig.  9.  (a)  m  vs.  pressure  for  various  blade  inlet  setting  angle  and  (b)  m  vs.  efficiency  for  various  blade  inlet  setting  angle. 


x  104 


Fig.  10.  m  vs.  losses  for  various  blade  inlet  setting  angle. 


Fig.  11.  Comparison  between  the  VTB  simulation  results  and  experiment. 


results  presented  in  the  reference  paper  [16].  The  dimensions 
used  for  the  model  are  same  as  that  of  the  experimental  compres¬ 
sor.  As  shown  in  the  figure,  the  VTB  simulation  results  match 
very  well  with  experiment  in  the  neighborhood  of  the  surge  line. 
At  the  same  time,  the  deviation  begins  to  increase  as  mass  flow 
rate  goes  up.  The  reason  for  this  might  be  due  to  some  other 
losses  (disc  friction,  choking,  etc.),  which  are  not  taken  account 
in  our  model. 

Figs.  12  and  13  show  the  comparison  between  the  VTB  sim¬ 
ulation  results  and  the  corresponding  results  presented  in  the 
reference  paper  [1].  The  agreement  between  the  results  is  very 
good,  with  VTB  slightly  over  predicting  the  efficiency  and  the 
pressure. 


Table  3 

Centrifugal  compressor  parameters  used  for  comparing  with  the  experimental 
results 


Parameter 

Value 

Parameter 

Value 

Blade  inlet  angle  (°) 

40 

Inlet  hub  diameter  (m) 

0.054 

Vane  inlet  angle  (°) 

28 

Outlet  tip  diameter  (m) 

0.18 

Inlet  tip  diameter  (m) 

0.106 

Outlet  hub  diameter  (m) 

0.18 

Fig.  12.  Comparison  between  VTB  simulation  results  and  paper  [1]  results  in 
terms  of  efficiency. 
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Fig.  13.  Comparison  between  VTB  simulation  results  and  paper  [1]  results  in 
terms  of  pressure. 


Table  4 


Parameters  for  dynamic  centrifugal  compressor  model 


Parameter 

Value 

Parameter 

Value 

Slip  coefficient 

0.84 

Outlet  tip  diameter  (m) 

0.06 

Diffuser  inlet  angle  (°) 

70 

Outlet  hub  diameter  (m) 

0.06 

Inlet  tip  diameter  (m) 

0.0465 

Mean  impeller  channel  length  (m) 

0.025 

Inlet  hub  diameter  (m) 

0.0158 

Mean  diffuser  channel  length  (m) 

0.006 

output  of  this  system  is  electrical  power.  The  heat  power  output 
dissipated  into  the  air  is  not  considered  here.  The  outlet  pres¬ 
sure  is  kept  at  120,000  Pa  in  the  simulation.  It  can  be  pointed 
out  her  that  the  compressor  efficiency  remains  almost  constant 
because  the  change  in  mass  flow  rate  is  very  small.  The  simula¬ 
tion  parameters  for  the  compressor  are  shown  in  Table  4. 

6.  Conclusion 

In  this  paper,  an  analytical  model  for  the  centrifugal  compres¬ 
sor  is  developed  from  the  first  principles  where  energy  transfer  is 
being  taken  into  consideration.  The  model  developed  allows  the 
user  to  predict  the  compressor  performance  from  the  geometric 
information.  In  this  paper,  the  incident  loss  and  friction  loss  are 
modeled,  so  are  other  losses:  clearance  loss,  backward  loss  and 
volute  loss.  A  dynamic  model  for  the  compressor  is  then  pro¬ 
grammed  into  the  VTB  simulation  environment  as  a  component 
of  an  electrical  system.  From  that,  we  can  predict  the  compres¬ 
sor  performance  curves  such  as  outlet  pressure,  efficiency  and 
losses.  In  addition,  surge  line  obtained  from  the  simulation  result 
can  be  used  to  define  stable  operation  range.  The  dynamic  sim¬ 
ulation  model  can  be  used  as  a  virtual  test  bed  for  compressors. 

Appendix  A 


5.4.  System  efficiency 

Fig.  14  shows  total  system  efficiency  versus  mass  flow  rate 
and  compressor  efficiency  versus  mass  flow  rate.  The  effective 


Mass  flow  rate  (kg/s)  Xl  o3 


Fig.  14.  System  efficiency. 


A.l.  Velocity  triangles  for  flow  through  the  stage  of  a 
centrifugal  compressor 


Velocity  diagrams  for  flow  through  the  stage  of  a  centrifugal 
compressor  are  shown  in  Fig.  A.l. 

Fig.  A. la  refers  to  the  inlet  velocity  triangle  and  (b)  to  the 
outlet  velocity  triangle.  All  the  variables  defined  in  both  (a)  and 
(b)  will  be  used  for  the  following  energy  transfer  and  losses 
analysis. 

The  absolute  inlet  velocity  C\  is  given  by 
1 

Ci  =  - m 

PoiM 


where  poi  is  the  constant  stagnation  inlet  density  and  Ai  is  the 
inlet  reference  area.  The  tangential  velocity  of  the  impeller  U \ 
is  given  by 


D  i 
2 


.  co  =  7i  D\ N 


where  N  is  the  rotational  speed  of  the  impeller  and  D\  is  the 
average  diameter,  which  is  defined  as 


1 

2 


Ai  +  dIo 
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not  have  a  tangential  component  of  velocity. 


A/lideal  =  —r  =  ^02  •  U2  (A.3) 

m 

The  ratio  between  Cq 2  and  U2  in  a  radial  vane  impeller  is  called 
the  slip  factor,  which  is  defined  as 

C02 

0  =  - 

U2 

From  above  analysis,  the  ideal  specific  enthalpy  delivered  to  the 
fluid  is 


Wc  o 

A/Jideal  =  -=a-U\  (A.4) 

m 


U2 


where  D{  1  is  the  inducer  tip  diameter  and  Dhi  is  the  hub  casing 
diameter. 

A.  2.  Energy  transfer  for  an  ideal  compressor 

For  the  ideal  compressor,  the  ideal  specific  enthalpy  delivered 
to  the  fluid  will  be  derived  as  follows. 

The  torque  equals  the  change  in  angular  momentum  of  the 
fluid  in  an  ideal  compressor, 

rc  =  m  •  (r2Co2  fiCe  1)  (A.l) 

The  power  delivered  to  the  fluids  is 

Wc  =  cd  •  rc  =  co  •  m  •  (r2Ce2  -  nCe  1) 

=  fn  •  (U2Cq2  -  U\Cq\)  =  m  •  A/qdeal  (A.2) 

In  order  to  simplify  the  analysis,  we  assumed  that  there  were  no 
stationary  pre- whirl  vanes  and  air  approaching  the  impeller  did 
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